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We consider the tunneling of localized excitations (many boson bound states) in the presence 
of a bosonic bath. We show both analytically and numerically that the bath influence results in 
ry^ • a dramatical enhancement of the amplitude of the excitation tunneling. The order of the bosonic 

flow in the course of the tunneling process is obtained. On the background of the giant tunneling 
enhancement we observe and describe additional resonant enhancement and suppression of tunneling 
due to avoided level crossings. 
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O ■ I. INTRODUCTION 



Tunneling properties of, e.g., local bond excitations in molecules, of bound states in quantum lattices or of large 
spins in the presence of anisotropy fields share a lot of similarities of corresponding model approaches and technical 
treatments. Despite the amount of published work this field lacks complete understanding. Especially intriguing are 
the relation to so-called dynamical tunneling in phase space [|l[ , the influence of heat bath coupling g,^ , the exact 
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■ O ■ computation of tunneling probabilities [Qj^] , to name a few. 

Q ' To study different aspects of the problem it is helpful to formulate a model which shares most of the needed 

o ; . . . 

properties, and is numerically solvable (because it is typically hard to obtain exact analytical results). 

^~~^ In the present paper we will study the tunneling of multi boson bound states in a model which we believe posesses 

0^ ' all necessary qualitative ingredients needed for the understanding of a more general case. 

^^ , The Hamiltonian of the model is given by 

a^■ 

^^ ' 1 r 2 21 

O' -^ — o ("]^'^i) +io-to.2) + C (a^ a2 + aia2) + S (a^ as + aia^ + a2 as + a2a^) (1.1) 

O' 2^ ^ 

jrt ' Here the operators ai,ai {I = 1,2,3) are bosonic annihilation and creation operators with standard commutation 



relations [a;, a^] = Sim- Note that (1.1) is invariant under permutation of site 1 with site 2. 



y , Equation (El) describes, e.g., a three site (or bond) molecule, i.e. a trimer [j. It could also serve as a model for 



O ' tunneling in a dimer [p| (just sites 1 and 2) under the simultaneous influence of a heat (or energy) bath on site 3, 

^ ] coupled to the symmetric modes of the dimer. It may be used for a study of excitations in H2O molecules, where the 

S^ ■ two O — H stretching modes (dimer) are additionally coupled to the relative angle between the two bonds (third site) 

Finally (f]^) may be transformed into a model of a spin with local anisotropy and external field (site 3) coupled to 
a bosonic field (using Schwinger bosons for sites 1 and 2). The coupling is such that the total spin length will not be 
conserved, but may fluctuate. The dimer limit 6^0 may be represented as a spin length conserving model with the 
Hamiltonian H, = S^, + 2CS^ 10. 



Model (1.1) conserves in addition to the energy also the total particle number b which is an eigenvalue of the total 



number operator B 

B — ni+n2 + n^ , ni—afai . (1.2) 

Note that due to this conservation the infinite dimensional Hilbert space decomposes into an infinite set of finite 
dimensional orthogonal subspaces, each with a given total particle number. Consequently numerical evaluations will 
always translate into diagonalizations of finite matrices, without worrying about truncation errors. The classical limit 
is obtained by replacing the operators ai with complex scalars Vt; (Hermitian conjugation is replaced by complex 
conjugation). 

Besides the above mentioned general interest in such kinds of model we want to give some additional reasons 



for studing of (LI). During the past decade we witnessed an explosion of interest to a new class of excitations 
in lattice models. These models are typically described by Hamiltonian equations of motion of degrees of freedom 
associated with lattice points of spatial periodic lattices. In the presence of interaction between these degrees of 
freedom, when considering linear equations of motion, the translational symmetry of the underlying lattice allows 
one to obtain delocalized plane wave solutions (phonons, magnons and other whateveions. The spectra of these 
excitations have finite upper bounds due to the discreteness of the system. If instead considering rather arbitrary 
additional nonlinearities in the equations of motion, a new generic type of solutions appears: discrete breathers 
pT[ . These discrete breather solutions are time-periodic spatially localized excitations and appear in one-parameter 
families. The discrete breather concept can be easily extended to dissipative systems |l2j. 

Much more complicated is the quantization of these excitations. The question arises which eigenstates of the 
quantum Hamiltonian operator correspond to classical discrete breathers. Intuition leads to the conclusion that one 
should search for multi whateveron bound states [13[ . One signature should be the appearance of extremely narrow 
quasiparticle bands in the high energy sector. 'Quasiparticle' stands for the need of just one wave vector which labels 
the states of the quantum breather band (as opposed to many wave vectors which label many whateveron continua) 
Q,|5[. This view is supported by the consideration of a (semiclassical) tunneling process of a discrete breather, 
which implicitely assumes that the quantum system may support coherent tunneling of localized excitations similar 
to a (quasi)particle |l^. 'Extremely narrow' could have different meanings. In the best case these bands should be 
much narrower than the mean level spacing in the given energy sector. This seems to be not so simple owing to the 
increase of the density of states when going up in energy. Another aspect of 'narrow' means weak interaction with 
nearby lying states. 

In contrast to the classical case, where implicit and explicit existence proofs of discrete breathers are known [[l7|-p0[ , 
and where reliable computational tools are available wich allows one to obtain numerical solutions on large lattices 
pi[, the situation with the quantum case is different. Except for some integrable one-dimensional models no analytical 
solutions are known. Numerical evaluations will have to face the tremendous problem of diagonalizing large matrices 
even for moderate system sizes (like 6-10 sites). In order to gain an insight into the problem it is useful to study 



small systems where the discrete breather concept can be partially used. Although there is no localization length 
anymore, we may still find excitations which are mainly located on one site. The smallest possible system is thus a 



two site system, i.e. a dimer. For (5 = equation (1.1) describes an integrable dimer. The analysis of this model M 
demonstrates that the existence of classical discrete breather solutions corresponds to the appearance of tunneling 
pairs of eigenstates. These pairs of eigenstates have energy splittings which are orders of magnitude smaller than the 
averaged level spacing. The relative splittings decrease with the increase of the eigenenergies, in full correspondence 
with the above described expectations. 

One nongeneric feature of the dimer model is its integrability. Although it allowed to obtain analytical results, it 
could hide some features of the generic nonintegrable case. This conclusion led to the consideration of the trimer 
model. The third site was deliberately chosen to be different from the first two sites in order to preserve the symmetries 
of the model. Some classical and quantum properties of ( |l.l| ) were described in detail in a previous paper ||7| . The 
discussion of the classical properties was devoted mainly to periodic orbit bifurcations which allow for an appearance 
of discrete breather solutions transforming into each other under the 1 ^^ 2 permutation in the trimer. In the quantum 
case we again observe pairs of symmetric and antisymmetric states with nearly coinciding energy levels. Note that the 
energy splitting in such pairs is determined by the phase space tunneling amplitude between the above two classical 
solutions. 

The paper m] presents data on the behavior of such pairs as well as of individual levels with varying coupling 
parameter S. We observe that pairs can be crossed by other pairs and/or by individual levels. The crossing for levels 
with the same symmetry is certainly avoided whereas levels with differing symmetries really cross each other. In many 
cases a pair is able to survive many crossings and restore itself up to large values of S. However, sooner or later a 
fatal crossing comes which destroys the pair. In the classical system this loss is manifested by a strong overlap of the 
corresponding quantum wave function with chaotic layers in the phase space of the classical model. 

The level splittings were not analized in more detail. Instead only their behaviour on scales comparable with the 
averaged level spacing were of interest. Here we will present a careful analysis of the mighty changes of the tiny 
splittings as the coupling S is changed. Denote with x,y,z the eigenvalues of the site number operators ni,ri2,n.3. 
We may consider the quantum states of the trimer at (5 = when z is a good quantum number and then follow the 
evolution of these states with increasing S. The state ioi S — can be traced back to C = and be thus characterized 
in addition by x and y. 

II. TUNNELING ORDER, MAGNITUDE AND RESONANCES 

In Fig. 2(b) in the previous paper |[7| we picked a tunneling pair which has energy w 342 for 5^0. The fixed 
model parameters were B = 40 and C = 2. The pair states are characterized by a; = 26(0), y ~ 0(26) and z = 14 for 
C = S = 0. Note that this pair survives approximately 30 avoided crossings before it finally is destroyed at coupling 
strength S w 2.67. 

Figure 1 presents the dependence of the tunneling splitting, AE, for this pair. One sees that the splitting rapidly 



increases gaining about eight orders of magnitude when S changes from to sUghtly above 0.5. Then this rapid but 
nevertheless smooth rise is interrupted by very sharp spikes when the sphtting AE rises by several orders of magnitude 
with S changing by mere percents and then abruptly changes in the opposite direction sometimes even overshooting its 
prespike value. Such spikes, some larger, some smaller, repeat with increasing 6 untill the splitting value approaches 
the mean level spacing (of order one in the figure) . Only then one may say that the pair is destroyed since it can be 
hardly distinguished among the other trimer levels. 
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FIG. 1. Level splitting versus 5 for level pair as described in the text. Solid line - numerical result. Dashed line - semiclassical 
approximation as described in text. Filled circles - location of wave function analysis in Fig. 2. 



Another observation is presented in Figs. 2 (a-e). We plot the intensity distribution of the logarithm of the squared 
symmetric wavefunction of our chosen pair for five different values oi 6 — , 0.3 , 0.636 , 1.0 , 1.8 (their locations 
are indicated by filled circles in Fig.l). We use the eigenstates of B as basis states. They can be formally represented 
as \x, y, z > where x, y, z are the particle numbers on sites 1, 2, 3, respectively. Due to the commutation of B with H 
two site occupation numbers are enough if the total particle number is fixed. Thus the final encoding of states (for 
a given value of 6) can be chosen as |a;, z) (see also [Q for details). The abscissa in Figs. 2 is x and the ordinate is 
z. Thus the intensity plots provide us with information about the order of particle flow in the course the tunneling 
process. 



Fig. 2a, Flach et al 



Fig.2b, Flach et al 




Fig. 2c, Flach et al 






Fig.2d, Flach et al 
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Fig.2e, Flach et al 




FIG. 2. Countour plot of the logarithm of the symmetric eigenstate of the chosen tunneling pair (cf. Fig.l) for five different 
values oi 5 — 0, 0.3, 0.636, 1.0, 1.8 (their location is indicated by filled circles in Fig.l). 
In Fig. 2(a) three equidistant grid lines are used. In Figs.2(b-e) ten grid lines are used. Minimum value of squared wave function 



is 10" 



maximum value is about 1. 



For 6 — (Fig. 2a) the only possibility for the 26 particles on site 1 is to directly tunnel to site 2. Site 3 is decoupled 
with its 14 particles not participating in the process. The squared wavefunction takes the form of a compact rim in 
the (x, z) plane which is parallel to the x axis. Nonzero values of the wavefunction are observed only on the rim. This 
direct tmmcling has been quantitatively described in pi. When switching on some nonzero coupling to the third site, 
the particle number on the dimer (sites 1,2) is not conserved anymore. The third site serves as a particle reservoir 
which is able either to collect particles from or supply particles to the dimer. This coupling will allow for nonzero 
values of the wavefunction away from the rim. But most importantly, it will change the shape of the rim. We observe 
that the rim is bended down to smaller z values with increasing 6. That implies that the order of tunneling (when, e.g., 
going from large to small x values) is as follows: first, some particles tunnel from site 1 to site 2 and simultaneously 
from site 3 to site 2 (Fig. 3(a)). Afterwards particles flow from site 1 to both sites 2 and 3 (Fig. 3(b)). With increasing 
d the structure of the wavefunction intensity becomes more and more complex, possibly revealing information about 
the classical phase space flow structure. 
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FIG. 3. Order of tunneling in the trimer. Filled large circles - sites 1 and 2, filled small circle - site 3. Arrows indicate 
direction of transfer of particles. 



Thus we observe three intriguing features. First, the tunneling splitting increases by eight orders of magnitude 
when b increases from zero to 0.5. This seems to be unexpected, since at those values perturbation theory in 6 should 
be applicable (at least Fig. 2b in |j demonstrates that this should be true for the levels themselfes). The second 
observation is that the tunneling begins with a flow of particles from the bath (site 3) directly to the empty site 
which is to be filled (with simultaneous flow from the filled dimer site to the empty one). At the end of the tunneling 
process the initially filled dimer site is giving particles back to the bath site. Again this is an unexpected result, 
since it implies that the particle number on the dimer is increasing during the tunneling, which seems to decrease 
the tunneling probability, according to the results for an isolated dimer. These first two results are closely connected, 
as we will show below. The third result concernes the resonant structure on top of the smooth variation in Fig.l. 
The resonant enhancements and suppressions of tunneling are related to avoided crossings. Their presence implies 
that a fine tuning of the system parameters may strongly suppress or enhance tunneling which may be useful for 
spectroscopic devices. In the following we will explain all three observations. 

A. The order of tunneling 

To understand the order of tunneling we consider a generalized model Hamiltonian for C = 5 = 

i^g = /("!) + 7(^2) +5('^3) . (2.1) 



It reduces to (1.1) when choosing ]{x) ~ x^ /2 and g{z) = 0. Switching on some interaction inside the dimer (C ^ 0) 
we will observe some tunneling probability. The smallness of this probability is due to the presence of an energy 
barrier hindering successive transferring of particles from one site of the dimer to the other. This energy barrier is 
given by 

Ab= max |/(5-zo)-(/(6-zo-x) + /(x))| (2.2) 

0<2;<(6-zo) 



where the total number of particles b and the number of particles on the still isolated third site zq are fixed. Note 
that the top of the barrier is reached when the particles are evenly spread inside the dimer (it is straightforward to 
generalize to other cases). 

At nonzero S the value of z may fluctuate. Now we have to consider the energy in the two-dimensional plane (x, z). 
The energy landscape for an arbitrary distribution of particles with respect to the chosen initial distribution x — 0, 
y = b — Zq, z = Zq rcads 

Afa(x, z) = \fix) + f{b-x-z) + g{z) - ,f{b - zo) - g{zo)\ . (2.3) 

For small values of 6 the variation of z will be small and we may expand around zq: z = zq + Sz. As noted above 
( p. 31) has a maximum at a; = (6 — zq)/2 for fixed zq. In general the derivative of this function along z will, however, 
not vanish in this point. This implies that a variation of z will change the barrier height. We have to determine the 
sign of the variation 6z which lowers the barrier height. This sign will tell us whether the particle number on site 3 
will decrease or increase during the tunneling. Thus, we may predict or explain the change of the rim shape of the 
above described squared wavefunction. This in turn will allow us to conclude about the order of tunneling. Define 
ei — f{b — Zq) — 2/((6 — zo)/2) and 62 = f'{{b — zo)/2) — g'{zo). Straightforward arithmetics leads to the result: 

5z > if [ei > and 62 < 0] or [ei < and 62 > 0] (2.4) 

Sz < otherwise (2.5) 



Our chosen Hamiltonian (l.I) yields Sz < which is in full accord with the observed change of the rim shape in Figs. 
2. 

The above considerations allow us to conclude that the tunneling enhancement due to the lowering of the barrier 
is not sensitive to the sign of Sz. The only possibility to avoid tmmcling enhancement is to choose 62 = 0. Still one 
will be faced with consideration of higher order contributions of Sz to the barrier change. 

B. Quantitative account for the tunneling enhancement 

Now we present a semiclassical approach which is capable of accounting for the rapid increase of the energy splitting 
in the tunneling pair shown in Fig.l. For this sake we shall first discuss tunneling in the dimer. Then we will show 
how the introduction of an interaction with the third site will modify the results. 

A semiclassical treatment of tunneling in the dimer was outlined in Jl^l where canonical coordinates were introduced, 
using the classical variables ^/ — y^e*'''' , 

q=-[n2-ni], p ^ (p2 - fi, 

J = ni+n2, X = 5 (<*5i + '/52) • (2.6) 

These coordinates allow one to represent the dimer Hamiltonian in the form 



^d=U) +'^' + 2CW('0 -q^cosp (2.7) 

Here the phase p plays the role of a momentum whereas q is the coordinate along which the tunneling in the dimer 
takes place. The coordinate J is an integral of motion, since the Hamiltonian does not depend on its canonically 
conjugated momentum x- Iii order to calculate the tunneling probability amplitude we first choose two symmetry 
breaking isolated periodic orbits with a certain energy E and then look for a classical path ^(t) in the complex time 
plane which connects these two orbits. On this way we have to cross classically forbidden areas where the momentum 
p becomes imaginary. What we need to know is the dependence of the momentum p on the coordinate q which, in 
the classically forbidden area, is determined by the equation 



E^[^) +ql + 2cJ (^] - ql coshp. (2.8) 



The boundaries ±qo{E) of this area or the turning points arc obtained from equation (2.8) at p = for a given energy 
E. 

The tunneling probability amplitude can be estimated as 

W ~ expj-S*} 

where 

S= J p{q)dq. (2.9) 

-90 



This integral can be readily evaluated for the energy E^ — J^/2 + C^ (for a given value of J). Then go = \J J"^ 1^ — C^ 
and 

2C 



S{E = E,) = -^J^^AC^ + J\n{j—j===^ (2.10) 

This expression coincides, to within the preexponential factor, with the results obtained in |lq| . However our aim 
here is to consider the tunneling amplitude of the trimer or, to be more specific, consider how the tunneling splitting 
for a certain pair of the trimer levels varies with the increase of the coupling parameter 5. 

A naive approach is based on equation ( |2.10| ) . An interaction with the third site of the trimer inay lead to a change 



of the number J of particles in the dimer (first two sites). One can readily see that the action ( 2.1C ) becomes larger 
for larger J values. Since the total number of particles in the trimer is fixed, & = J + 713, an increase of J occurs 
only at the expense of the third site which may provide the necessary particles. Therefore the tunneling amplitude 
decreases if n-^ should decrease, and this is completely opposite to what we have observed numerically in the previous 
section. It means that a more detailed analysis is needed. 

First, we introduce canonical coordinates which arc explicitly symmetric with respect to permutation of the first 
two sites of the trimer. 



^ -{ni+n2)~ -113, TT = ^{(pi + (P2) - (P3, 



b = ni+n2 + n3, x = M'/'i + </'2 + V's)- 



(2.11) 



The definitions of tlie two coordinates q and p, related to tlie dimer, remain as in equation (2.6). Then the trinier 



Hamiltonian (1.1) takes the form 



H= {}b+ iQ) +q^ + 2Cd{\b + \Q) ~ q^ cosp 



26M\b - Q) {lb +\Q + q) cos {it + \p) 



2Sj{\b - Q) (lb +^z-q) cos (tt - ip) 



(2.12) 



It becomes straightforward to obtain equations of motions for the four relevant coordinates, with b being an integral 
of motion and x redundant. We need here only two of this equations 



Q= - 25J{\b - Q) {\b + iQ + g) sin {tt + \p) 



28J{\b-Q){\b+\Q-q)s\n{'K~\p) 



(2.13) 



ig 



C{\b+lQ) 
{U+\Qf~q^ 



: cosp - 



5 cos (tt 



5 cos (tt — ip) 



Ib+^Q + q) 



{¥-q) 



:\b-Q) 



'\b+\Q-q) 



(2.14) 



These equations will allow us to estimate the deviation of the tunneling trajectory in the (g, Q) plane. The latter 
can be directly mapped onto the (z, x) plane. Again we may consider trajectories connecting isolated periodic orbits 
with p = p = TT = 'K = Q. A reasonable approximation is to assume that tt = tt = along the path and then use 



equation ( 2.1g ) to estimate the deviation of the Q value from its initial Qo values, 

.^ A At VQo(6/3-Qo) 

^^^""^ 2 ^ ^C 

Here Qmax is the largest value of Q. It is estimated by substituting Q = Qq, q = (the middle point of the tmmcling 



(2.15) 



path) and Pmax into equation (2.13). As for Pmax it is found for the given dimer energy from equation (2.7) at g = 



At ~ l/-\/J2/4 — C^ is the time of motion through the forbidden region (traversal time). 
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First, we again obtain that the tunnehng trajectory really deviates towards lower values of Q, i.e. towards a 
decreasing number z of bosons in the third site. It fits our numerical results. In order to calculate the corresponding 
action we slightly deform the trajectory. Better to say we straighten the trajectory, which finally simply becomes a 
straight line with constant z which is shifted towards lower z-values. We believe that this deformation would lead to 
a negligible change of the result since the deformation is carried out mainly in the borders of the forbidden area. The 
principal contribution to the integral comes, however, from the central part of this area. 

This procedure allows us to return to the calculation of the tunneling amplitude in the dimer by means of the 



integral (2.9). However, the number of bosons, which has been J in the isolated dimer, becomes now J' = J + AQ. 



As for the energy of the dimer it should remain without change {E — E^), since we neglect small corrections (order 
of (5^) due to the interaction with third site. The distance to the upper bound £^3 ~ J'^/2 + C* can be readily found, 
AE « JAQ. There should be also a small shift of the turning points q'q which is, however, of a minor importance. 
As a result we can represent the action in an integral form 

"" :^ + AE-q^ 

dq arccos ^ , (2.16) 

rj72 ^ ^ 

.0 ^C\l--<i' 



This integral can be calculated analytically and in the trivial case of AQ = it reduces to equation ( ^.10 ). However 



the general analytic expression for (ref3-ll) is extremely cumbersome. That is why we have preferred to calculate this 



integral numerically. The thick dashed line in Fig.l is the graphical form of ( 2.16 ). We observe very good agreement 
with the numerical data for tunneling in the trimer up to 6 values where possibly our expansion becomes too crude 
and where the spiky structure sets in. Yet our consideration is capable of quantitatively describing the raising of the 
tunneling splitting by eight orders of magnitude! 

C. Resonances due to avoided crossings 

Let us give a semiquantitative understanding of the spiky behaviour of the energy splitting dependence on S in 
Fig.l. Comparing the position of the spikes with the level variation in Fig. 2b in H we immediately realize that each 
resonance in Fig. 1 is related to an avoided crossing of our level pair with either a single third level or another tunneling 
pair. A simple analytical description of such crossing can be obtained following the same simple reasoning which has 
led to the understanding of an avoided crossing of just two levels (see, e.g., pl[|). 

Let us first assume that we have found all eigenenergies and eigenfunctions of the trimer for given values of the 
parameters C and 6q. We denote the energies and wave function of our tmmcling pair as eis, £ia and ^pis, tpia- Here 
indices s and a relate to the symmetric and antisymmetric states. We also assume that there is a level or another 
pair of levels which happen to be very close to our tunneling pair. Then we may consider small variations S of the 
coupling parameter S near its Sq value (5 = 6q + S) and calculate the variation of the level energies by means of the 
perturbation theory, neglecting contributions from all other more distant levels. The corresponding perturbation is 
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explicitly symmetric with respect to the 1 — 2 permutations. Therefore it may mix only the states with the same 
symmetry, so that those levels will avoid crossings. As for states of different symmetries, their crossing is allowed. 

1. Pair — single level crossing 

First we consider the case when the tunneling pair is approached by a single level which we choose to be symmetric 
with respect to the permutation of the dimer sites. The energy and the wave function of this level a,i 5 = 5q are £2s 
and '(/'2s- Such a crossing can be described by the 3x3 matrix 

( eis + Vsii~5 V,i2~5 \ 

Vsi2~5 e2s + K22 . (2.17) 

V eia + Vall) 

where all terms proportional to S are the corresponding matrix elements of the perturbation. They vanish at Sq which 
is close enough to the avoided crossing. 

Diagonalizing the matrix ( p. 17 ) we obtain three levels 



^^1,2} = 7) ± V -. + l/sl2'^' 



Ea = eal + Vall6 (2.18) 

There are two different scenaria which depend on the parameters of the matrix. In the first case at S = Sq the 
order of the levels in either increasing or decreasing order is Sia, ei^, £25- In this case we will first see a decrease 
of the splitting down to zero (!) and a subsequent increase of the splitting up to some maximum with subsequent 
decrease down to its original value. In the second case £2^, £la, £ls the above sequence is reversed - we first observe 
an increase up to a maximum, then a decrease down to zero and a final increase to the original splitting value. All 
the parameters may be chosen to be positive. The difference between the coefficients Vgn and Vaii is neglected since 
this difference reflects the slow change of the level splitting with varying 5 (cf. Fig.l). The perturbation theory allows 
us to deal only with small values of S, hence parameters should satisfy the conditions that \esi — £s2| ^ \Vsii — Vs22\ 
and |14i2| < 1. 

Then the tunneling splitting A can be found as a function of d. It can be defined as the smallest distance between 
the antisymmetric level and the two symmetric levels, 

A = min{|^,i - Ea\, \E,2 ~ Ea\} (2.19) 
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FIG. 4. Level splitting variation at avoided crossings. 
Inset: Variation of individual eigenvalues participating in the avoided crossing. Solid lines - symmetric eigenstates, dashed lines 
- antisymmetric eigenstates. 

An example is shown in Fig.4(a) (T^n = Kii = 1,V;22 = -3, T4i2 = 0.001, e,i == 0.1, e,2 = 0.3, eai-£.i = 0.00002, 
the asymmetric level lies belovif its symmetric counterpart in the tunncUng pair.) where A is plotted as a function of 
I. 

Crossing of a tunneling pair with an antisymmetric level is described in a similar manner by the same matrix with 
interchange of indices s and a. 

2. Pair — pair crossing 

In order to analyze a pair — pair crossing we should consider two tunneling pairs with the energies Sgi , £ai, and 
£s2i £a2 together with the corresponding wave functions, ipsi, ipai, and "082, '0a2- Again we assume that these two 
pairs lie close to each other at (5 = (5o, so that the perturbation theory can be applied, disregarding all the other states 
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of the trimer. Then we should consider the matrix 

V^sl2^ £2s + Vs22 





V 

whose diagonahzation results in four energies 



\ 



Val25 ea2 + ^22^ / 



(2.20) 



esl+£s2 + {Vsll+Vs22)~5 ^ 
^s{l,2} - ^ ± 



'[esi-es2 + {Vsii-Vs22)5\ 



VA26' 



(2.21) 



^a{l,2} 



£al + £a2 + [Vail + Va22)5 ^ , / [g gl - £a2 + (Kll - ^22)^]^ , ,,3 ?2 

^ + ''012'' 



± 



We may now construct four possible differences between the energies of two symmetric and two antisymmetric states 
which give us possible splittings. Then choosing the difference with the lowest absolute value we define the tunneling 



splitting A at a given value of the parameter 5. The choice of other parameters in equations (2.21) is governed by 
the similar principles as in the pair — single level crossing. 

There are at least three different scenaria which may appear. The first two are realized when the interaction between 
the two symmetric and the two antisymmetric states is much stronger than the interaction inside each pair which 
yields the individual pair splittings. If Vsi2 ~ Vai2 the avoided crossing results in a simple resonance of the splitting 
(Fig. 4(b) , parameters £,i ^ 0.1, £,2 - 0.3, £„! = 0.09999, £a2 = 0.2999, Vsn = Van = 1, ^^s22 = K22 = -3, 
^sl2 = Val2 = 0.001). 

If the symmetric states are interacting stronger (or weaker) than their antisymmetric counterparts Vsi2 7^ ^112, 
two intersections inside one of the two pairs are possible (Fig. 4(c), parameters £^1 — 0.1, £s2 = 0.3, £ai = 0.09999, 

£a2 = 0.2999, V^ii = Vail = 1, K,22 = ^22 - -3, ^.12 = 0.002, Kl2 = 0.001). 

Finally, if the interaction between the two symmetric and the two antisymmetric states is much weaker than the 
interaction inside each pair, the pairs will intersect each other and the avoided crossing structure will appear on a 
finer scale than the pair splittings. An example is shown in Fig. 4(d) (parameters £^1 = 0.1, £^2 — 0.3, £qi — 0.09, 

£a2 = 0.2, y,n = Vail = 1, 1^22 = ^22 = -3, ^^12 = Val2 = 0.1). 



III. CONCLUSION 

We have described and explained peculiar features of dynamical tunneling of a strong excitation in a trimer. 
The tunneling amplitude may be enhanced due to the presence of the third site with fiuctuating particle number. 
Numerical calculations show that this amplitude grows about eight orders of magnitude when the coupling parameter 
6 still remains much smaller than all other parameters in the system. An interesting feature of this process is that the 
third site serves as a donor which injects bosons into the dimer and thus speeds up tunneling. The injected bosons are 
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of course returned to the third site at the end of the tunnehng process. Such a process thus needs a nonzero amount 
of bosons to be deposited on the third site at the beginning of tunnehng. 

It is "worthwhile to compare the observed enhancement of tunnehng with our consideration of a tunnehng of a 
breather in a one-dimensional chain mm . The result of that paper was that the interaction of the dimer with other 
degrees of freedom in the chain may result only in a reduction of the tunneling amplitude. This finding was in 
agreement with studies of disspative tunneling PJ3| . The reason for the opposite result in the present work is that we 
consider a situation in which the third harmonic site is excited to a rather high level and is able to provide bosons 
necessary to enhance the tunneling. This is in contrast to the above references where the temperature was assumed 
to be zero or low, meaning that the bath degrees of freedom were not excited. 

Note that the discussion in Subsection II. A prompts a possibility of a tunneling enhancement even at zero occupation 
of the third site. This can be achieved in our model by choosing 3(^3) — 0^3713 with lu^ > B/2. Again the expected 
enhancement of tunneling is opposite to the above statement and is due to the choice of such a large value of the 
frequency W3 of the third site which makes the adiabatic approximation ||l6| invalid. It means that the shape of the 
instanton calculated in |l6| may be significantly altered by the interaction with this high frequency mode. 

The spikes observed in the S dependence of the tunneling amplitude are well explained by the relevant crossings 
with cither single levels or with other tunneling pairs of levels. At higher values of 6 when the tunneling splitting 
comes close to average level spacing (1 in our model) the whole picture becomes more complicated and one may think 
about an increasing role of chaos assited tunneling pi. However, discussion of this process goes beyond the framework 
of this paper. 
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